Transport properties of partially equilibrated quantum wires 
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We study the effect of thermal equilibration on the transport properties of a weakly interacting 
one-dimensional electron system. Although equilibration is severely suppressed due to phase-space 
restrictions and conservation laws, it can lead to intriguing signatures in partially equilibrated 
quantum wires. We consider an ideal homogeneous quantum wire. We find a finite temperature 
correction to the quantized conductance, which for a short wire scales with its length, but saturates 
to a length-independent value once the wire becomes exponentially long. We also discuss thermo- 
electric properties of long quantum wires. We show that the uniform quantum wire is a perfect 
thermoelectric refrigerator, approaching Carnot efficiency with increasing wire length. 

PACS numbers: 71.10.Pm, 73.23.-b 



o 
o 



> 

On 
OO 

m 
m 

O 

&\ 

o 



X 



I. INTRODUCTION 

The quantization of the dc conductance in ballistic 
quantum wires, first observed about two decades ago^ 1 ^ 
is one of the fundamental discoveries of mesoscopic 
physics. The staircase-like dependence of the conduc- 
tance on the electron density, with plateaus at integral 
numbers of 2e 2 /h is readily understood from the single- 
electron pictured The latter associates each plateau with 
a fixed number of occupied electronic subbands, each sup- 
plying one quantum of conductance 2e /h. On the other 
hand, interactions between one-dimensional electrons of- 
ten lead to qualitatively new phenomena. These are 
commonly described within the so-called Luttinger liq- 
uid theory, 4 drastically different from Landau's Fermi liq- 
uid description applicable to higher-dimensional systems. 
The remarkable success of the simple single-electron pic- 
ture in describing the quantization of conductance is at- 
tributed to the fact that quantum wires are always con- 
nected to two-dimensional leads, where interactions be- 
tween electrons do notplay a significant role. In fact, it 
was shown in Refs. that in an ideal Luttinger liq- 

uid connected to Fermi liquid leads, the dc conductance 
is completely controlled by the latter and, therefore, is 
not affected by interactions in the wire. 

For that reason, the discovery of small 
temperature-dependent deviations from perfect 
quantization§i2ii2iiiii2ii2iiliSii& of the conductance 
of quantum wires at low electron densities raised a lot 
of interest. These generally manifest themselves as a 
shoulder-like structure just below the first plateau of 
conductance. Weak at the lowest temperatures available, 
this feature becomes more significant as the temperature 
is increased, turning into a quasi-plateau at about 
0.7 x (2e 2 /h). A number of theoretical efforts trying to 
reveal the microscopic mechanism of this so-called "0.7 
structure" have been made. Several spin-related ap- 
proaches attribute the effect to spontaneous polarization 
of the electron spins in the wir8 & 17 i 18 or the existence 
of a local spin-degenerate quasi-bound state playing the 



role of a Kondo impurity. 9 ' u Other approaches discuss 
the role of scattering from plasmonsj^i spin waves^ or 
phonons^ 

Despite the absence of a commonly accepted micro- 
scopic theory, it is generally recognized that electron- 
electron interactions must be included to account for 
the effect. As a consequence, a number of recent 
publications reconsider the effect of interactions on 
the transport properties of one-dimensional conduc- 
tors, going beyond the picture of an ideal Luttinger 

liquid3»25 l 26 l 27 1 28,29 1 30 1 3l l 32 1 33,34,35 l 36 1 37 Here we focus Qn 

a very fundamental aspect of interactions, studying how 
they lead to the equilibration inside the wire of electrons 
coming from the two leads. We emphasize that this effect 
is absent in an ideal Luttinger liquid. Indeed, the bosonic 
elementary excitations of the Luttinger liquid have infi- 
nite lifetime, thus there is no relaxation towards equilib- 
rium in these systems, no matter how strong the interac- 
tions. Within the Luttinger- liquid theory the processes 
leading to the equilibration of the electron system would 
be accounted for by the additional terms in the Hamil- 
tonian, which are irrelevant in the renormalization group 
sense. Instead of pursuing this strategy, we consider the 
regime of weakly interacting electrons, thereby avoiding 
the complexity of the Luttinger-liquid picture. 

Non-interacting electrons propagate ballistically 
through the wire and, therefore, keep memory of the 
lead they originated from. Thus the distribution func- 
tion of electrons inside the wire depends on the direction 
of motion. For the right- and left-moving particles it is 
controlled, respectively, by the left and right lead: 
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Here e p is the energy of an electron with momentum p and 
6{p) is the unit step function. The left and right leads are 
assumed to have the same temperature T, but different 
chemical potentials /i; = fi + eV and fi r — fi (see Fig. [TJ. 
Using the distribution function ([T]) one easily finds the 
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FIG. 1: Schematic picture of the quantum wire of length L 
which is formed by confining a two-dimensional electron gas 
with gates (dark regions) . Electrons in the left and right lead 
are described by Fermi distribution functions characterized 
by temperature T and chemical potentials fii and fi r , respec- 
tively. 
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which coincides with the well-known conductance quan- 
tum 2e 2 /h up to an exponentially small correction 

In the presence of interactions, the ballistic propaga- 
tion of electrons through the wire may be interrupted by 
collisions with other electrons. As a result of these col- 
lisions, some electrons change their direction of motion 
thus losing the memory of the lead they originated from. 
Such scattering processes modify the electron distribu- 
tion function which is then no longer given by Eq. |T|). 
The effect of the electron-electron collisions on the dis- 
tribution function depends strongly on the length of the 
wire. Indeed, electrons traverse short wires relatively 
fast, so the interactions do not have the time to change 
distribution dXJ) considerably. On the other hand, in the 
limit of a very long wire one should expect full equili- 
bration of left- and right-moving electrons into a single 
distribution, even in the case of weak interactions. 

To simplify the subsequent discussion, in this paper we 
consider the case of electrons with quadratic spectrum, 
€ p = p 2 /2m, where m is the electron effective mass. Then 
the system is Galilean invariant, and one can easily in- 
fer the electron distribution function in the fully equili- 
brated state. Viewed from a frame moving with the drift 
velocity Vd — I/ne (where I is the electric current and 
n is the electron density) the electron system is at rest 
and must be described by the equilibrium Fermi distri- 
bution. Performing a Galilean transformation back into 
the stationary frame of reference this distribution takes 
the form 



fp 



1 
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(3) 



where the chemical potential \x and temperature T in- 
side the equilibrated wire are, in general, different from 
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FIG. 2: (a) Energy-conserving two-particle scattering process 
violates conservation of momentum. Such processes can occur 
only in inhomogeneous systems, (b) Dominant three-particle 
collision which gives rise to corrections to the conductance 
of short quantum wiresi^ A hole at the bottom of the band 
scatters off electron excitations close to the Fermi level. 



fiur and T. At zero temperature, T = T = 0, the distri- 
butions ([1]) and ([3]) coincide, provided fj,u r = ft ± VdPF, 
where pf = Tthn/2 is the Fermi momentum of the sys- 
tem. At non-zero temperature the distribution function 
© of electrons inside the wire is slightly different from 
the distribution (QJ supplied by the leads. In a previous 
worl*2& we have shown that the mismatch of the distri- 
bution functions inside a very long wire and in the leads 
results in additional contact resistance, reducing the con- 
ductance to 
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It is worth noting that the quadratic in T / /x correction in 
Eq. ^ is much more significant than the exponentially 
small correction in Eq. 

The mechanism of equilibration of the electron distri- 
bution function in one dimension is not fully understood. 
While in higher dimensional systems equilibration at low 
temperature is primarily provided by pair collisions of 
electrons, these do not provide a relaxation mechanism 
in one dimension. This is due to the conservation laws 
for momentum and energy which severely restrict the 
phase-space available for scattering processes, Fig. [U[a). 
As a result, pair collisions in one-dimensional wires can 
only occur with a zero momentum exchange or an in- 
terchange of the two momenta, leaving the distribution 
function unaffected. The leading equilibration mecha- 
nism thus involves collisions of more than two particles. 
For a weakly interacting system, it is then natural to 
assume that equilibration is provided by three-electron 
scattering processes. 

The effect of three-particle collisions on the transport 
properties of short wires has been studied in a recent 
work by Lunde, Flensberg, and Glazmani22 In such short 
systems the effect of equilibration is weak and the dis- 
tribution function can be calculated perturbatively from 
the distribution of non-interacting electrons |1| within 
the Boltzmann equation framework. Following this ap- 
proach, Lunde et al£^ obtained interaction-induced cor- 
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rections to transport, which they attributed to specific 
three-particle scattering events that change the number 
of left- and right-movers. Indeed, in the absence of inter- 
actions, the current flowing through the system can be 
viewed as the superposition of the right- and left-moving 
flows of electrons supplied by the left and right leads, 
respectively. Once interactions are included, these indi- 
vidual contributions change due to electron-electron colli- 
sions, and one needs to account for the fact that electrons 
can now change direction. The electric current flowing 
through the wire is thus given by the sum of the non- 
interacting part Iq = GqV , and the change in, say, the 
number of right-moving electrons inside the wire 

I = G V + eN R . (5) 

Interaction-induced corrections to transport therefore 
arise from processes which change the number of right- 
and left-moving electrons rather than a change in the 
velocity of the charge carriers, as also pointed out in 
Ref. [23. 

As shown by Lunde et al.^ the most efficient three- 
particle process changing the number of right-moving 
electrons involves scattering of an electron into an empty 
state near the bottom of the band, see Fig. [2Kb) . By cal- 
culating the resulting N R , they obtained the correction 
to the conductance ^ of the wire of the form 

5G = - 2 -f±e-^, (6) 

where the length l eee is determined by the interaction 
strength and shows a power-law temperature depen- 
dence. The exponential smallness of the correction (J6j) 
is due to the small probability of finding an empty state 
near the bottom of the band. Since the small backscat- 
tering probability grows linearly with the length L, the 
correction SG oc L. 

Because the same three-particle processes are respon- 
sible for the thermal equilibration of the distribution ([1]) 
into © , the papers Ref. [26| and Ref. [13] reviewed above 
study the same physical phenomenon, albeit in the op- 
posite limits of a long and a short wire. In the present 
paper we bridge the gap between these two limits. We 
discuss how the electron distribution evolves from the 



out-of-equilibrium form (fTJ) in a short wire to a fully equi- 
librated form {3| in a long wire, and study how transport 
is affected by the process of equilibration. Our analysis 
focuses on weak electron-electron interactions. It is thus 
formulated entirely in terms of electrons, and does not 
use the bosonization technique. 

The paper is organized as follows. In sections HTIaiid Hlll 
we investigate how the conductance changes with increas- 
ing length of the wire. In Sec.[TT]we expand on the kinetic- 
equation treatment^! of backscattering in short wires and 
study the length dependence of the conductance while 
the correction SG remains exponentially small. In sec- 
tion lllll we turn to the regime of exponentially long wires, 
where the correction 5G ~ (e 2 //i)(T/^) 2 , cf. Eq. (j4|). In 
section HVl we study the thermoelectric effects and show 
that the uniform quantum wire is a perfect thermoelectric 
refrigerator, attaining Carnot efficiency with increasing 
wire length. Details of some calculations can be found in 
the Appendices. 



II. CONDUCTANCE OF SHORT WIRES 

Consider a quantum wire of length L, connected by 
ideal reflectionless contacts to non-interacting leads bi- 
ased by voltage V. We are interested in the process of 
thermal equilibration of the electrons inside the wire, i.e., 
in how the transition from distribution (fTJ to (|3|) occurs, 
and how it affects the transport properties of the system. 

Following Lunde et al.^* we describe the electron 
transport in the wire in the framework of the Boltzmann 
equation 

POh* =IpM ( 7 ) 
m ox 

We consider the steady-state setup in which the elec- 
tron distribution function f PiX depends on the position 
x along the wire, but not on time. The collision inte- 
gral Ip X [f] is, in general, a nonlinear functional of the 
distribution function, whose form is determined by the 
scattering processes affecting the distribution function. 
As discussed above, in our case the dominant processes 
are three-particle collisions, in which case 



IprAf) = - E E ^123;l'2'3'[/l/ 2 / 3 (l-/l0(l-/20(l-/3')-/l'/2'/3'(l-/l)(l-/2)(l-/3)], (8) 

C2 ,0-3 (Tj^crj, CT3 



where iui23;i'2'3' is the rate for scattering the set of in- 
coming states {picri,P20'2,P^s} into the set of outgoing 
states {p' 1 <y i ,p 2 <7 2j?'3 cr 3}i an d for notational convenience 
we shortened fa = f PilX - 



The Boltzmann equation ([7]) should be solved with the 
boundary conditions stating that the distributions f p0 of 
the right-moving electrons (p > 0) at the left end of the 
wire and f p .L of the left-moving electrons (p < 0) at the 
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right end coincide with the distribution function fp ^ in 
the leads, Eq. |T]). The conductance of the wire can then 
be found from Eq. ([5]), with the rate of change in the 
number of right-moving electrons related to the collision 
integral via 



N R = 



2 io dx r ^ ipAf] ' 



(9) 



Solving the Boltzmann equation exactly is a very dif- 
ficult problem due to the non-linearity of the collision 
integral l[B]). so one generally has to make some simplify- 
ing assumptions. Such assumption in our case is that the 
temperature T is small compared to the chemical poten- 
tial /i. 

Clearly, at T = no real scattering processes are al- 
lowed, and the unperturbed distribution f p — 9{pp — \p — 
mu\) solves the Boltzmann equation ([7]) for any value of 
the drift velocity u. Since in this case the collision inte- 
gral Ip, x [f] vanishes, we get N R — 0, and, according to 
(O, the conductance of the wire is 2e 2 /h. 

A finite temperature T acts in two important ways. 
First, it affects states near the Fermi level: the step in 
the zero-T distribution softens, providing partially oc- 
cupied states in a momentum range bp ~ T/vp around 
the Fermi points. Secondly, it ensures a finite occupa- 
tion of a hole (i.e., a vacant state) near the bottom of 
the band. Although the occupation probability of such 



a hole is exponentially small, 1 — f p 



, its pres- 



ence is crucial for the three-particle processes that change 
the number of right-moving electrons, see Fig. H^b) . It 
is important to realize that the backscattering of holes 
is accompanied by scattering of electrons near the Fermi 
points, Fig. [5Jb). In fact, this is the mechanism of the 
equilibration of the distribution function to the form 
in long wires. Although the backscattering rate is expo- 
nentially small, N R oc e~ M / T , it scales with the length of 
the wire. Thus the full equilibration is achieved in wires 
whose length L exceeds an exponentially long equilibra- 
tion length Z cq oc e M / T . The exact definition of l cq will be 
given below, see Eq. ([59")) . 

In this section we will discuss the case of short wires, 
L <^l eq . The regime L > Z eq will be discussed in Sec. IIIII 



A. Very short wires 

We start our discussion with the case of very short 
wires, recently considered by Lunde, Flensberg, and 
Glazmani^l The authors argued that for short enough 
wires, the interactions have little time to change the dis- 
tribution function from its initial value fp ^ given by 
Eq. |[T]). allowing one to perform a perturbative expan- 
sion in the scattering rate iui23;i'2'3'- In the lowest order, 
this amounts to approximating the collision integral as 



Ip,x[f] — Ip,x[f^ 



(10) 



function inside the wire, which they used to compute the 
electric current to first order in the scattering rate. 

The resulting correction to the conductance of the wire 
has the form © , in which microscopic details of the inter- 
action potential are absorbed into the length l eee . Lunde 
et aii22 performed their calculation for a specific model of 
electrons interacting via a potential defined by its Fourier 
transform V q = Vq(1 — q 2 /q^). This expression results 
from the expansion of a general potential under the as- 
sumption that small- momentum scattering is dominant. 
The parameter qo <^ kp accounts for the screening by the 
nearby metallic gates, while Vq is the zero-momentum 
Fourier component of the screened Coulomb potential. 
Within this model, the length l eee is given bj*2£ 



r 1 

( < ( 



V k F 



9o ) \M> 



(11) 



Solving the Boltzmann equation to this approximation, 
they obtained an expression for the modified distribution 



A more careful treatment of the Coulomb interaction 
screened by a gate leads to an additional logarithmic tem- 
perature dependence in Eq. (TlTj) . see Appendix |A"1 

To better understand the result ([6|) and find the limits 
of its applicability, we discuss the qualitative picture of 
this phenomenon. Let us focus on a single three-electron 
collision process. The most favorable collision involves 
a maximal number of states close to the Fermi points. 
However, due to the conservation of both energy and mo- 
mentum, collisions that change the number of right- and 
left-movers cannot occur near the Fermi level, and have 
to involve states deep in the electron band. As pointed 
out by Lunde et al.2L the scattering process most sus- 
ceptible to alter the current thus typically scatters two 
electrons close to the Fermi points and one electron at 
the bottom of the band, as schematically depicted in 
Fig. Hlb). It is convenient to think of this collision as a 
process in which a deep hole, corresponding to the outgo- 
ing electron state, is backscattered by electron excitations 
close to the Fermi level. These excitations are typically 
associated with a momentum change \5p\ ~ T/vp due 
to Fermi-blocking, so that the backscattering occurs over 
a distance ~ T/vp in momentum space. Let us further- 
more characterize this process by introducing a scattering 
rate 1/to, which can be approximated by a constant since 
the initial and final states both lie at the bottom of the 
band. 

The change N R in the number of right-moving elec- 
trons per unit time, due to these three-particle collisions 
can then be readily obtained. It is given by the product 
of the scattering rate 1/to for one such collision times 
the number of deep holes susceptible to be backscattered. 
The latter can be estimated from the probability to find 
a left- or right-moving hole / T and the number 

of states (T/vf)/(Ii/L) available within the typical mo- 
mentum range of the backscattering process. Taking into 
account that the scattering of a left- or a right-moving 
hole both contribute to N R , but with a different sign, 
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one finally has 



N R = l( e -» R / T -e-» L / T ) 



TL 

hvF 



.l^EL e -»/T L 

To flVF 



(12) 



where A/i 



fi L , and we absorbed potential nu- 



merical prefactors into the definition of to. Throughout 
this paper we use subscripts I and r to denote the left 
and right leads, whereas superscripts L and R refer to 
the left- and right-moving electrons. In short wires the 
chemical potentials of electrons are not significantly af- 
fected by the scattering processes, so fi R — /i;, fi L = fx r , 
and A/i = eV. 

We then notice that according to Eq. ([5]) the correction 
to the conductance of the wire due to the backscattering 
processes is SG = eN R /V. As a result we recover the 
result ^ of Lunde et ai.jSl provided 



'cee 

v F ' 



(13) 



The derivation of Eq. ^ relied on the assumption that 
the occupation probability of a deep hole is well described 
by the distribution of non-interacting particles, or, alter- 
natively, that one can approximate the collision integral 
according to Eq. (|10[) . This approximation holds in cases 
where the hole typically scatters no more than once dur- 
ing its propagation through the wire, and any transition 
between subsystems of left- and right-movers occurs in a 
single collision. One thus expects this result to be valid 
for wires shorter than the mean free path of the hole Iq. 
Since the typical momentum of a hole contributing to N R 
is of order T/vp, we estimate Iq ~ Ttq/pf- Substituting 
the estimate ([IB"]), we obtain 



o 



(14) 



For the particular model of the interaction potential used 
in Ref. 123 we estimate 



Vokp \ fkp^ 



(15) 



In wires longer than Iq, holes near the bottom of the 
band experience multiple collisions while they propagate 
through the wire, the distribution function deviate sig- 
nificantly from the unperturbed form ([1]), and the result 
© is no longer applicable. 



B. Longer wires: lo <C L <C l e q 

In wires longer than a typical hole near the bottom 
of the band is scattered many times while traversing the 
wire. Each collision changes its momentum by a small 
amount Sp ~ T/vp <C pf, with a sign that varies in a 



random fashion. The hole thus performs a random walk 
in momentum space. This picture is analogous to the 
diffusion of a Brownian particle in air. In the latter case, 
the change of momentum of the particle in each collision 
is small because its mass is much larger than that of the 
air molecules. Similarly to the case of Brownian motion, 
one can use the small parameter Sp/pp ~ T/fj, to bring 
the collision integral of holes to a much simpler Fokker- 
Planck form 



d ( Id 
Ip,x[9l ~~qZ \ A(p)9p,x ~ ^-Q-[B{p)g p , x ] 



(16) 



where we introduced the hole distribution g p x = 1 — f p x . 
The functions A(p) and B(p) entering Eq. (fT6|) are model 
specific. In the case of three-electron collisions they can 
be determined explicitly. They depend on the three- 
particle scattering rate as well as the electron distribution 
function in the vicinity of the Fermi level. The latter can 
be assumed to be unperturbed by the collisions in the 
wire as long as L <C Z cq oc e^l T . The resulting derivation 
of A(p) and B(jp) can be found in Appendix |Bl here we 
provide order of magnitude estimates. 

First we notice that B(p) has the physical meaning 
of the diffusion coefficient in momentum space, i.e., the 
typical momentum change of a hole over time t behaves as 
(Ap) 2 ~ Bt. Assuming as before that the hole changes its 
momentum by ±T/vp once during time ro, we conclude 
that (Ap) 2 ~ (T/vp) 2 t/T for t > t . Thus we estimate 



B 



T 2 



VpTQ V F l e 



(17) 



where we used our earlier estimate (|13[) of t , and the 
microscopic expression for l eee is given by Eq. (|11[) . 

Although B is a function of momentum p, the typical 
scale of the variations of B(p) ispF. Thus for the particle 
at the bottom of the band one can approximate B(p) by 
its value at p — 0, which we will denote as B. Then A(p) 
is easily obtained by noticing that the collision integral 
(fT6| has to vanish if the hole distribution function takes 
an equilibrium Boltzmann form 



„(0) _ p p 2 /2mT -m/t 



(18) 



This condition leads to the relation A(p) = Bp/2mT, 
which is also confirmed explicitly in Appendix [Bl Using 
this result one easily transforms the Boltzmann equation 
to the form 



V dg P ,x 
m dx 



B d 
~2dp 



mT 



9p,: 



dp 



(19) 



The boundary conditions express the fact that the dis- 
tributions of the right-moving holes at the left end of the 
wire and that of left-moving holes at the right end are 
controlled by the respective leads: 



e p 2 /2mT e -(p+eV)/T 



e p 2 /2mT e -ix/T 



for p > 0, (20a) 
for p < 0. (20b) 
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Here we again assumed \i r = fi and fii = fx + eV. The 
kinetic equation in the form (|19[) is applicable only to 
exponentially rare holes with \p\ <C Pf- Thus the Fermi 
statistics of the holes is irrelevant, and the boundary con- 
ditions on the distribution function have the Boltzmann 
form. Finally, combining our earlier results ([5|), ((9]), and 
(|16|) . we express the correction to conductance of the wire 
as 



with 



N R 



B 
J 



dx 



d 9. 



p,x 



dp 



(21a) 



(21b) 



p=0 



i.e., conductance is determined by the behavior of the 
distribution function near p = 0. 

The solution of equation (fTi?)) with boundary condi- 
tions ([20)) shows two different regimes, depending on the 
length of the wire. In relatively short wires the effect of 
hole scattering is weak, and to first approximation one 
can assume that the distribution function g PtX does not 
depend on position x and coincides with the distribution 
([2H)l provided by the leads. This distribution is discon- 
tinuous at p = 0, namely g p — > e _Al! r / T at p — > ±0. 
To be more precise, one should notice that the Fokker- 
Planck approximation applies to wires of length in the 
range Iq <C L <C l cq . At the lower end of this range, 
L ~ lo the holes near the bottom of the band are scat- 
tered a few times by the electrons near the Fermi level 
and change their momentum by ~ T/vp. Thus in the 
center of the wire the discontinuity of the distribution 
g p is smeared by (Ap)o ~ T/vp. As the wire length 
increases, the diffusion of holes in momentum space be- 
comes more pronounced, and at a certain length scale 1% 
the smearing Ap reaches a larger scale (Ap)i — (mT) 1 ' 2 . 
(Indeed, (Ap) 1 /(Ap) ~ y/Jt/T > 1.) We shall consider 
the regimes L <C ^i and L l\ separately, as different 
approximations can be applied to the kinetic equation 
(fT9f in these two cases. The estimate for the length scale 
l\ will be obtained below, see Eq. (|24|) . 



1. Wires of length in the range < L < li 



Let us present the hole distribution function as g 



g(0) 

and 



g, where g PtX is the equilibrium distribution (|18j) 
l]X is the correction caused by the applied bias V. 

(At small bias we expect g oc V.) The distribution g p °l, 
of course, satisfies the kinetic equation (fl9|) . Then, since 
equation (| 19jl is linear in g, it also fully applies to g P , x . 
It is important to note, however, that the two terms in 
the right-hand side of this equation are not of the same 
order of magnitude. Indeed, at p ~ Ap we have dg/dp ~ 
g/Ap > (p/mT)g, provided Ap < (toT) 1/2 . Thus the 
propagation of holes through the wires of length L in the 
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FIG. 3: As one goes along the wire from the left lead to the 
right one, the distribution function of holes changes. The 
number of left-moving holes decreases as some of them turn 
around and start moving to the right. 



range Iq <C L <C h is described by the simplified equation 

P dg P ,x _ B d 2 g P:X 



m dx 



2 dp 2 



(22) 



To find the correction to conductance (|2ip for a wire in 
the regime ^ « L « I[ one needs to solve this equation 
with the appropriate boundary conditions deduced from 
Eq. (|2"0|) . We leave such a complete solution for future 
work. Instead, we perform a simple dimensional analysis 
to conclude that the step in the distribution function near 
p = is broadened by 



Ap ~ (BmL) 



1/3 



(23) 



This result can also be obtained from a simple physical 
argument. Figure [3] shows the hole distribution function 
g p ^ x at different positions along the wire. The scatter- 
ing processes contributing to the electric current involve 
holes entering from the right lead with momentum Ap. 
moving to the left with their velocity gradually decreas- 
ing, and eventually returning to the right lead. In order 
to lose the momentum of order Ap the hole has to ex- 
perience sufficiently many collisions in the wire, which 
requires time t determined from the standard diffusion 
condition (Ap) 2 ~ Bt. Propagating through the wire at 
a typical velocity Ap/m until the turning point, the hole 
will move by distance (Ap/m)t ~ L. Combining these 
two estimates, we recover our earlier result (|23[) . 

At this point we can estimate the upper limit l\ 
on the length of the wire L, to which the approach 
used here is applicable. In order to neglect the term 
B(d/dp)[(p/2mt)g] in the right hand side of equation 
we assumed Ap < (raT) 1/2 . From Eq. we see 
that this approximation fails when the length of the wire 
reaches the value 



h 



(mT 3 ) 1 / 2 
B 



1/2 



— I 

T 



(24) 



where we also used our earlier estimate (fTTf of B. As 
expected, l\ 3> Iq, see Eq. (TT4"|) . i.e., the approach used 
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in this section applies to a parametrically broad range of 
wire lengths. 

To find the effect of three-particle scattering on the 
conductance of the wire we use the boundary conditions 
(l20l) to estimate 



dg p . x 
dp 



p=0 



Ap T 



Substituting this estimate into Eq. (|2"Tj) , we obtain the 
correction to the conductance in the form 



2/3 



-/i/T 



(25) 



This correction should be compared with the result ([6]) of 
Lundc et al. 27 Both expressions are exponentially small 
and grow with the length of the wire, but correction ([23)1 
shows a slower growth, 8G oc L 2 / 3 , rather than linear 
growth in Eq. ©. One can easily check that at the 
crossover, L ~ Iq, the results and (|25p are of the 
same order of magnitude. 



2. Wires of length in the range h -C L -C ?eq 

As mentioned above, a hole near the bottom of the 
band performs a random walk in momentum space. In 
the case of wires longer than l\ one needs to carefully 
consider the effect of the parabolic spectrum of the hole 
—p 2 /2m. This spectrum plays the role of a potential 
barrier for the random walker, see Fig. [4] In order for 
the hole to backscatter, and thus change sign of its mo- 
mentum, it has to overcome the barrier. The rate of 
such backscattering events is controlled by the height of 
the barrier measured from the Fermi level, and is ex- 
ponentially small as e~ p / T . Evaluating the prefactor of 
backscattering rate is an interesting problem, similar to 
that of a Brownian particle escaping from a local mini- 
mum of the external potential. The general features of 
this problem are well understood^ In order to overcome 
the barrier the particle has to not just reach the top, but 
move beyond it far enough for the potential to drop below 
the maximum by more than the temperature T. Applied 
to a hole diffusing in momentum space, this means that 
the bac ksca ttering is controlled by the region of width 
Ap ~ V mT around p = 0. 

In wires shorter than l\ this process cannot fully de- 
velop because of the small time needed for the hole to tra- 
verse the wire, and one has to include into consideration 
the spatial dependence of the distribution function g P}X . 
At L ^> l-y the holes spend enough time inside the wire 
to fully complete the backscattering process. Therefore, 
away from the ends of the wire the distribution function 
no longer depends on position. As a result, the left-hand 
side of the kinetic equation (fP9|) vanishes, and it can be 
rewritten in the form 



d_ ( P_ | dg p 

dp \ mT ® p dp 



= 0. 



(26) 












P 


"/'/ / 




\ /'/ 







V(x) 



FIG. 4: A hole with spectrum — p 2 /2m performs a ran- 
dom walk in momentum space in the presence of a barrier. 
Backscattering of a hole is analogous to a Brownian particle 
escaping from a local minimum of the potential, see inset. 



To complete the mathematical formulation of the prob- 
lem one has to impose the appropriate boundary condi- 
tions on the distribution function. Since equation (|26[) ig- 
nores the spatial dependence of the distribution function, 
we cannot reuse the boundary conditions (|2"D|) . Instead 
we assume that the chemical potentials of the right- and 
left moving holes are established by the leads and do not 
vary along the wire: 



eP 2 /2mT e -( M +eV)/T for p ^ 
9P - i eP V2mT e - M /T 5 for _ p>> ^ i 27 i 



This assumption implies that the total number of 
backscattered holes is too small to affect the chemical 
potentials. This is justified by the fact that the backscat- 
tering rate is exponentially small. In wires of exponen- 
tially large length L > Z eq this condition is violated. The 
latter regime will be discussed in Sec. IIIII 

The solution of equation (|26|) with boundary condi- 
tions (12~T|) is straightforward and gives 



e p72mT e - M /T 



1 



:- eV / T - 1 
v27rmT 



z-v' 2 ' 2mT dp' ) . (28) 



This distribution function smoothly interpolates between 
the boundary conditions (|27[) imposed by the applied 
bias. As expected, the crossover occurs in a narrow re- 
gion of width Ap ~ \JmT at the bottom of the band. 
To linear order in eV/T we find 



d9p 



eV 



-PIT 



V2irmT 3 
resulting in the backscattering rate 

eVBL 



N R = 



ti/T 



h\ // 2TrmT 3 



(29) 



(30) 



see Eq. (|21b|) . As a result, the correction to the conduc- 
tance (|21a[) takes the form 



5G 



2e^L 



-n/T 



(31) 
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where we have used the following precise definition 



h 



B 



(32) 



of the length l\, consistent with our earlier estimate ([24)1 . 
It is worth mentioning that for wires of length L ~ l\ the 
expressions ([23]) and (f3"Tj) give the same estimate for 8G. 

Our result (f3Tj) has the form similar to the prediction 
^ of Lunde et al£l for short wires, L <C lo. Both ex- 
pressions for the correction to the conductance are expo- 
nentially small, but grow linearly with the length of the 
wire L. However, due to the sublinear growth (|25|) in 
the intermediate range of wire lengths Iq -C L <C h, the 
prefactor l^ 1 in Eq. (f3"Tj) is parametricaUy smaller than 
l- e \ in Eq. ©, see 

The result (|3"Tj) can be derived qualitatively, following 
arguments similar to the ones used in Sec III Al There 
the change N R in the number of right-moving electrons 
per unit time was estimated as the ratio of the number 
of holes likely to backscatter and the average time To of 
such backscattering event, see Eq. (fT2"|) . Compared to the 
case of very short wires considered in Sec III Al for a hole 
to change direction, it must now cover a larger distance 
(Ap)i ~ V mT (Ap)o ~ T/vf in momentum space 
set by the smearing of the discontinuity of the distribu- 
tion function at the bottom of the band. The number of 
states available for the passage is thus larger by a factor 
(Ap)i/(Ap) ~ \J n/T compared to the case of a very 
short wire, Sec III Al On the other hand, even though the 
typical time between two three-particle collisions is still 
given by tq, it now takes many such collisions for a hole 
to complete the backscattering process. Because the hole 
performs a random walk in momentum space, the time 
Ti it takes to cover the longer distance (Ap)i can be 
estimated from t%/to ~ (Ap)^/(Ap)§ ~ ^/T. Combin- 
ing both effects we find that the correction to the con- 
ductance (f3Tj) should be smaller than ((6|) by a factor of 
^/n/T, in agreement with Eq. (f24l . 



III. CONDUCTANCE OF LONG WIRES 

In short quantum wires, the distribution function of 
electrons remains close to the unperturbed form {!]) pro- 
vided by the leads. The main change due to the processes 
of electron collisions occurs near the bottom of the band, 
with the discontinuity at p = being gradually smeared 
as the wire length L increases. Because the disconti- 
nuity affects only electrons deep below the Fermi level, 
the effect of collisions is exponentially small. In particu- 
lar, this enabled us to neglect the effect of collisions on 
the chemical potentials and assume that to first approx- 
imation the right- and left-moving electrons remain in 
equilibrium with the left and right leads, respectively. 

A much more significant change occurs in long wires, 
L ~ / q oc e^/ T , for which the exponential suppression of 
the equilibration effects is compensated by a large sys- 
tem size. Once the length of the wire becomes exponen- 



tially large, the relaxation of the electron system becomes 
significant, and eventually, at L ^> l eq the distribution 
function assumes the fully equilibrated form of Eq. ((3]). 
Unlike the relatively minor modification of the distribu- 
tion function in short wires, the difference between the 
distributions |T]) and (|3|) is not exponentially small and, 
more importantly, concentrated near the Fermi points, 
rather than at the bottom of the band. In this section 
we consider the conductance of the partially equilibrated 
wires, of length L ~ l cq . We start by discussing the form 
of the electron distribution in this regime. 



A. Electron distribution function in the case of 
partial equilibration 

Let us consider a segment of the wire, whose length AL 



is small compared to the equilibration length l cq cx e 



This condition implies that a typical electron with energy 
near the Fermi level passes through the segment without 
backscattering. On the other hand, AL is assumed to 
be sufficiently large for electrons to experience multiple 
three-particle collisions, which do not result in backscat- 
tering. Under these circumstances, the electron distribu- 
tion function in the segment will achieve a state of partial 
equilibration, in which the numbers N R and N L of the 
right- and left-moving electrons are not changed by col- 
lisions. The form of this distribution can be obtained 
from a general statistical mechanics argument. The mul- 
tiple collisions occurring in the system will maximize the 
entropy of the non-interacting electrons 



S = -2j2if P ^ /„ + (!- f P ) Hi - f P )] , 



(33) 



while preserving the total energy, momentum, N R , and 
N L , given by 



E = 


2 2^ e pfp' 
p 


(34a) 


P = 


p 


(34b) 


N R = 




(34c) 




p>0 




N L = 


2 $3 / P . 


(34d) 




p<0 





Subtracting from the functional (|33[) the expressions for 
the conserved quantities (|34ap ~ (|34d|) with the Lagrange 
multipliers 0, —f3u, —f3fi R , and —f3fj, L , respectively, and 
differentiating with respect to f p , we find that the max- 
imum of entropy is achieved for the distribution 



fp 



6{p) 



6(-p) 



3 (£„-«p-^ R )/T _|_ I e (e p -up-[t L )/T _|_]/ 



(35) 



Here T = 1/0 is the effective temperature, parameter u 
has dimension of velocity and accounts for conservation 
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of momentum in electron collisions, fi R and fi L are the 
chemical potentials of the right- and left-moving parti- 
cles. 

It is worth mentioning that the distribution (|35[) does 
not apply to particles near the bottom of the band, p ~ 
\JmT . Indeed, for a hole near p = collisions with and 
without backscattering (i.e., the change of the sign of p) 
are roughly equally likely. Thus the above discussion is 
not applicable in this case. In order to find the form of the 
distribution function near the bottom of the band, one 
should perform an analysis similar to that of Sec. Ill B 21 
In particular, the exponentially small discontinuity of the 
distribution (f35|) at p — will be smeared. On the other 
hand, most quantities of interest are determined by the 
behavior of the distribution function near the Fermi level. 
For instance, using (|35p we obtain the electric current in 
the form 

/ = + erut i (36) 

up to corrections small as e~~^/ T . Here Afi = /i R — fi L . 

It is instructive to see how the distribution (|35p inter- 
polates between the regimes of no equilibration |T]) and 
that of full equilibration ([3]). The unperturbed distribu- 
tion (fTJ) is obtained from (|35|) by setting u = and iden- 
tifying the chemical potentials with those in the leads: 
fi R = fii and ix L — /i r . In this case A/z = eV, and 
Eq. (|36|) reproduces the Landauer formula. The fully 
equilibrated distribution ((3| is obtained from (|35j) by set- 
ting Afi = [i R — fi L = 0. In this case the electric current 
(|36p is expressed as I = enu, which identifies parameter 
u with the drift velocity v$. 

In the regime when the distribution function (|35|) dif- 
fers from the limiting cases (|TJ) and ([3]) it is convenient 
to quantify the degree of equilibration in the wire by the 
parameter 

•» = -• (37) 

Vd 

The case of no equilibration corresponds to r/ — and 
that of full equilibration to r\ = 1. 

The meaning of the distribution function (|35|) can be 
further clarified by considering the Boltzmann equation 
([7]). The scattering processes contributing into the colli- 
sion integral ([5]) fall into two categories. The strongest 
processes preserve the numbers of the right- and left- 
moving electrons, whereas the ones resulting in backscat- 
tering are exponentially weak, as discussed by Lunde et 
al£L and also above in Sec. [TTJ Let us approximate the 
collision integral © by neglecting the weak backscatter- 
ing processes. Then, by substituting the distribution l|35p 
into the right-hand side of Eq. (jSJ), one easily sees that 
each term in the sum vanishes. Thus the distribution 
([55]) solves the Boltzmann equation ([7]) in this approx- 
imation. Furthermore, in the absence of backscattering 
the solution (|33|) applies for any choice of parameters T, 
u, fi R , and fi L , and in particular, for any degree of equi- 
libration rj. The value of r\ is ultimately determined by 



the exponentially weak backscattering processes and the 
length of the wire. 

B. Conservation laws 

Conductance of a long quantum wire, in which the elec- 
tron distribution function is fully equilibrated, was stud- 
ied in Ref. [26| , where a power-law correction to the quan- 
tized conductance was obtained, Eq. ([4]). The derivation 
of this result was based on an analysis of conservation 
laws for the number of electrons, energy, and momentum 
satisfied in electron collisions. Here we perform a similar 
analysis for a partially equilibrated wire. 

Conservation of the total number of particles N im- 
plies that in a steady state the particle current j(x) is 
uniform along the wire. Correspondingly, we infer from 
the conservation of total momentum P and total energy 
E that in the steady state a constant momentum current 
jp and a constant energy current je flow through the 
system. In the following it will be convenient to express 
these currents as the sum of the individual contributions 
from left- and right-moving electrons, e.g. j = j R + j L , 
thus introducing 



J R/L (x) = 


J — 00 


j^6(±p)v p fp^, 


(38a) 


j R P '\x) = 


J — 00 


-^e(±p)v p pf P! , 


(38b) 


iTV) = 


p oc 
J —00 


-?-6(±p)v p £ p f p/x . 


(38c) 



Here v p = p/m is the electron velocity, the positive sign 
in the step function corresponds to right-movers, while 
the negative one to left-movers. 

Near the ends of the wire, the distribution function of 
incoming electrons is controlled by the leads. Close to 
the left lead, the distribution f R of right-moving elec- 
trons thus assumes the form of the first term in Eq. |T|), 
and similarly, close to the right lead, the left-movers' dis- 
tribution f L is given by the second term in Eq. |1| . This 
allows us to readily calculate, for example, the current 
j R (l) of right-moving electrons near the left end of the 
wire. Unlike the total current j, the current j R (x) is not 
uniform throughout the system, since the equilibration 
processes ensure the conversion of right-moving electrons 
into left-moving ones. From conservation of the num- 
ber of particles it follows that the total number of right- 
moving electrons changing direction per unit time equals 
the difference between their outgoing and incoming flows 
at the right and left leads 

N R = j R (r) - j R (l). (39) 

A calculation of j R (r) requires the knowledge of the dis- 
tribution of right-moving electrons at the right end of the 
wire. As the latter is unknown, we proceed by express- 
ing j R (r) in terms of the total particle current j and the 
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incoming flow of left-movers supplied by the right lead 
3 R (r)=j-j L (r), (40) 

where j L (r) can now be determined from the known elec- 
tron distribution in the right lead. Combining Eqs. (1391) 
and P0"]) we can relate N R to the total incoming flow of 
particles as 

j R (l)+3 L (r)=j-N R . (41) 

Using the lead distribution function fl} the left-hand side 
of Pdj) is readily calculated, and takes the form GoV/e, 
with the conductance Go defined by Eq. ([2]). Noticing 
that the electric current / = ej, we then recover Eq. (|5j). 
For the purposes of this section we do not need to keep 
the exponentially small corrections to Go, and upon sub- 
stitution Go = 2e 2 /h we are left with the relation 

^V = j-N R (42) 

between voltage, current, and the rate of change of the 
number of right moving electrons due to collisions. 

Let us now analyze the consequences of energy con- 
servation in electron-electron collisions. Repeating the 
above steps for the energy current jg, we arrive at an 
expression 

im+3 L E {r)=j E -E R , (43) 

analogous to Eq. (|4"Tj) . Here E R is the rate of change of 
the energy of right-movers. 

The conservation of the number of electrons and energy 
ensure that the currents j and js are constant along the 
wire. It is convenient to combine them into the heat 
current 

3Q = 3E - W, (44) 

which is consequently also independent of position. 38 
Combining Eq. (03]) and (gT]) we find 

3§(!)+3Q(r)=jQ-Q R , (45) 

where 

Q R = E R — nN R (46) 

is the heat transferred into the right-moving subsystem 
by electron collisions. 

The left-hand side of Eq. (|45[) is the heat current in a 
non-interacting quantum wire. Direct calculation shows 
that it is exponentially small (see also the discussion in 
the beginning of Sec. IIV[) . and for our purposes here can 
be assumed to vanish. We therefore conclude 

JQ = Q R - (47) 

Since the heat current jq does not depend on position, 
it can be calculated at any point in the wire. In the 



regions not too close to the leads the distribution function 
is expected to have the partially equilibrated form (|35[) . 
Then, using the expressions (|38ap and (|38cp for j and j E , 
we obtain 

7T 2 T 2 

JQ = nu (48) 

6 \x 

to leading order in T I pi,. As expected, in the absence of 
equilibration, u = 0, the heat current vanishes. 

In Sec. IIII Al we introduced the distribution function 
(|35[) by discussing a short segment of the wire. The four 
parameters of this distribution T, u, /i R , and may, in 
principle, vary along the wire. The independence of heat 
current (|48|) on position then shows that the velocity u 
and, therefore, the degree of equilibration r\ are constant 
along the wire. Furthermore, since the electric current 
(|36[) and u are constant along the wire, one concludes 
that the difference of the chemical potentials A/j, = /i R — 
fi L is constant as well. The only two parameters of the 
distribution (|35|) that can vary along the wire are the 
temperature T and the average chemical potential (/j, r + 
ji L )/2. Their dependences on position are discussed in 
Appendix O 

C. Relation between the degree of equilibration r\ 
and the conductance of the wire 

To make further progress we elaborate on the relation- 
ship between the rates N R and E R , whose explicit forms 
depend on the details of the equilibration mechanism. As 
we discussed in Sec. IIII Al in the absence of scattering 
processes which change the number of left- and right- 
moving electrons, the distribution (|35|) is unaffected by 
electron-electron collisions, i.e., the collision integral ([5]) 
vanishes. In particular, for u = the unperturbed dis- 
tribution ([I]) supplied by the leads would retain its form 
in the wire. The backscattering processes, which by def- 
inition contribute to N R , also change the energy of the 
subsystem of right-moving electrons, resulting in a non- 
vanishing E R . Because both rates are caused by the same 
backscattering processes, one expects to find a relation 
between N R and E R . Here we establish such a relation 
with the help of conservation laws. An alternative and 
more formal derivation can be found in Appendix [Dl 

The backscattering processes transform the unper- 
turbed distribution fl]) to the partially equilibrated form 
(|35[) with non-vanishing u. The two distributions differ 
most prominently at energies within ~ T of the Fermi 
level. One can thus assume that all the right-moving elec- 
trons contributing to N R are removed from the vicinity 
of the right Fermi point and placed to the vicinity of the 
left one. Each such transfer reduces the momentum of 
the system by 2pp. Since the electron-electron collisions 
conserve momentum, a number of other electrons have 
to be scattered in the vicinities of the two Fermi points, 
see Fig. \5\ In the special case of three particle collisions, 
the transfer of electron from the right Fermi point to the 
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FIG. 5: Sequence of elementary scattering events leading to 
the transfer of one electron from the right to the left Fermi 
point. Conservation of momentum leads to the heating of the 
right-moving electrons and cooling of the left-moving ones. 



left one is accomplished in a number of small steps with 
momentum change Sp ~ T/vf, and at each step one ad- 
ditional electron is scattered near each of the two Fermi 
points, see Fig^b). 

As a result of the rearrangement of electrons near 
the two Fermi points, the momentum change 2pp of 
the backscattered electrons is distributed between the 
remaining right- and left-moving electrons, i.e., Ap R + 
Ap L = 2pp. Thus the energy of the remaining right- 
movers increases by AQ R = vpAp R whereas that of the 
left movers decreases, AQ L = —vpAp L . Then, the con- 
servation of energy requires Ap R = Ap L = pf- In the 
end, the energy balance for the right-moving electrons 
consists of a loss of ep due to removal of one particle 
from the Fermi level and a gain of AQ R — vppp = 2ep 
due to the redistribution of momentum. As a result, 
for every right-moving electron that changes direction, 
AN R = —1, the right-movers' energy increases by an 
amount AE R = ep ■ It is easy to check that the dif- 
ference between the chemical potential fi and the Fermi 
energy ep is irrelevant for our discussion, so we conclude 

-fiN R . 



E R 



(49) 



It is important to point out that this result is independent 
of a specific equilibration mechanism, or the degree to 
which equilibration has occurred. 

The result (|49|) can be also expressed in the form 

Q R = -2/ziV* (50) 

cf. Eq. ([4U)) , which expresses the simple fact that when 
a right-moving electron is moved to the left Fermi point, 
the remaining right-movers gain energy, see Fig.O Com- 
bining this expression with Eqs. (|47p and (|4"5)) . we obtain 



N R = -r nu, (51) 
where the dimensionless parameter ro is defined as 



Noticing that by definition of the degree of equilibration 
P7[) nu — rjnvd — rjj and using the conservation of the 
particle number in electron collisions, Eq. (|4"2")> . we find a 
linear relation between the applied bias and the electric 
current flowing through the wire. We can thus readily 
extract the expression for the conductance in partially 
equilibrated wires 



G- 2e2 



7T 2 T 



(53) 



where we discarded higher-order corrections in (T/ /i) 2 . 

The result (|53[) reaffirms that at finite temperature 
the processes of equilibration of the electron distribution 
function lead to a deviation of the conductance from its 
quantized value 2e 2 /h. In a fully equilibrated wire r] — ► 1, 
and this correction saturates at a value that does not de- 
pend on the details of the electron-electron interaction, 
reproducing our earlier result ((31) ■ The correction to the 
conductance in this situation is quadratic in temperature 
SG oc (T//i) 2 , in contrast with the results for the short 
wire, where SG ~ e~^l T . 

We obtained the expression (|53"|) for the conductance of 
the wire by taking advantage of the conservation laws for 
the electron-electron collisions as well as the basic prop- 
erties of the non-interacting leads the wire is connected 
to. Although the expression (|53p is thus very general, 
it does not fully determine the conductance of the wire, 
as the parameter r\ cannot be obtained in this approach, 
with the exception of special cases of non-interacting elec- 
trons, r\ = 0, and a very long wire, 7] = 1. To find r\ and, 
therefore, the conductance of the wire, for arbitrary wire 
length, one needs to consider a specific model of electron- 
electron interactions. We now turn to such a calculation 
for the most relevant case of relaxation via three-particle 
collisions. 



Partially equilibrated wires and equilibration 
length 



'■" =12 U 



(52) 



D. 



Our expression ([53")) for the conductance of partially 
equilibrated wires relied on the relation ([51]) between the 
rate N R of backscattering of right-moving electrons in 
the wire and the parameter u of the distribution func- 
tion ([3"5]) . Another relation between N R and the distribu- 
tion function can be found by considering the microscopic 
mechanism of such backscattering, using the approach of 
Sec. Ill B 2] Comparison of the two expressions will enable 
us to determine the degree of equilibration rj for wires of 
arbitrary length. 

The form of the electron distribution in a partially 
equilibrated state ([33)) is controlled by four parameters, 
namely, the temperature T, average chemical potential 
fj, = (/i R + [i L )/2, difference of the chemical potentials 
A/i = fi R — and the velocity u. In the absence of bias 
V applied to the wire, the temperature T and chemi- 
cal potential [i are equal to those in the leads, whereas 
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A/i and u vanish. As a result the distribution (|35|) re- 
produces the equilibrium Fermi-Dirac distribution, and 
clearly N R — 0. Our goal is to find N R to linear order 
in V. 

Applied bias affects the four parameters of the distri- 
bution (f35|) differently. As we discussed at the end of 
Sec. IIIIB| the temperature T and chemical potential /i 
acquire position dependence, whereas A/i and u no longer 
vanish, but remain constant along the wire. Thus, to lin- 
ear order in V one expects to find for the rate h R of 
backscattering per unit length of the wire 



L 



71 A/i + 7 2 u + 73 d x T + 74^/z. 



(54) 



Because the gradients d x T and d x fi are caused by bias 
applied to a long wire, they are not only proportional to 
V, but also scale as 1/L (see also Appendix |C|) . Thus 
for the exponentially long wires considered here, L ~ Z cq , 
effect of the gradients of T and /i can be neglected. It 
is also clear that 72 = 0. Indeed, at A/i = the distri- 
bution takes the fully equilibrated form ([3]) , for which no 
relaxation takes place, and N R = for any u. We thus 
conclude that the backscattering rate n R can be found 
for the simplest case of small A/i, vanishing u, and unper- 
turbed (position-independent) values of the temperature 
T = T and average chemical potential /_t. 

The resulting problem is equivalent to the one consid- 
ered in Sec. Ill B 21 The only difference is that because 
the length of the wire was assumed to be short, L <C l eq , 
the parameter A/i coincided with eV. Thus replacing 
eV -> A/i in Eq. flU) we find 



n r = -n 



2A/t 



Here the dimensionless parameter r\ is defined as 

L 



n = 



(55) 



(56) 



and the length l\ is given by Eq. (J32J). 

Our expression (|55p and the earlier result (f5Tj) relate 
N R to two different parameters of the distribution func- 
tion (|35|) . namely, A/i and u. Both of these parameters 
affect the electric current in the wire and can be expressed 
in terms of the drift velocity v c i = I/ne and the degree of 
equilibration 77. Indeed, according to the definition (|37|) 
of 77 and the expression (|36[) for the current, we have 



(57) 



2A/i 

nu = r}nvd, —r— = (1 - rjjnvd- 



Using these expressions to compare Eqs. (joTj) and ([55]) . 
we readily find the following expression for the parameter 

7} 

L 



V 



ra + n l cq + L 

where we have introduced the equilibration length 



L 



ro 



(58) 



(59) 



As expected, parameter 77 grows with the length of the 
wire from to 1. Because the backscattering process 
involves rare holes at the bottom of the band, the equi- 
libration length (jnij) at which the crossover occurs is ex- 
ponentially long at low temperatures. 

Substituting Eq. into ([53")) we find the following 
expression for the conductance of a quantum wire 



2e 2 



12 V>" 



(60) 



valid for L 3> l\. At L — > 00 it recovers the long wire 
limit while at l± <C L -C / cq it agrees with our earlier 
result J3T]). 



IV. THERMOELECTRIC PROPERTIES OF 
LONG WIRES 

We now turn to a situation in which the leads are not 
only biased by a finite voltage V but also exposed to 
a temperature drop AT. More specifically, we assume 
that the leads supply the wire with the following electron 
distribution 



(o) _ 



0(p) 



(-P) 



,(e p -»-eV)/(T+AT) + 1 e ( £p - M )/T + 1 ' 



(61) 



Our goal is to find the thermopower S, Peltier coefficient 
II, and the thermal conductance K of the quantum wire. 
These transport coefficients are defined by the following 
linear response relations 



V 

3Q 
3Q 



-SAT 



1=0' 



= ni 



AT=0' 



(62) 
(63) 
(64) 



The thermopower and Peltier coefficient are not indepen- 
dent properties of the system; they are connected by an 
Onsager relation II = ST. 

In the absence of interactions, electron distribution in 
the wire is given by Eq. (|6T|) , and its transport coefficients 
are easily understood. For instance, the thermopower 
and Peltier coefficient are given by 



n = st 



1 



■T(l- 



ln(l- 



,-m/t 



)]. (65) 



At low temperature T <C /i, the expression (|65|) is ex- 
ponentially small, II ~ (A*/ e ) e _M / T . The reason is that 
contributions to the heat current from electrons with en- 
ergies fj, + £ and /i — £ cancel each other, and the only 
reason II does not vanish completely is the absence of 
electronic states below the bottom of the band. In this 
paper we are not interested in such exponentially small 
results, unless the smallness can be compensated by a 
long length of the wire. Thus to first approximation the 
thermopower and Peltier coefficient of a non-interacting 
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quantum wire vanish. The latter conclusion can be easily 
obtained from the so-called Cutler-Mott formula 40 



S = 



ir 2 T d\uG 
3e dtp ' 



(66) 



generally applicable to systems of non-interacting elec- 
trons at T <§; tp- Considering that the conductance 
G = 2e 2 /h does not depend on the Fermi energy tp, 
one easily concludes that 5 = 0. 

To the same accuracy, i.e., neglecting corrections small 
as e _A1 / T , the thermal conductance of a non- interacting 
quantum wire is 



2tt 2 



(67) 



This expression can be derived by straightforward cal- 
culation of the heat current for the electron distribution 
(IBTj) with V = 0. Alternatively, one can obtain (f6"T|) from 
the Wiedemann-Franz law 



K 



3e 2 



TG, 



(68) 



by substituting the quantized conductance G = 2e 2 /h. 

It is important to note that both the Cutler-Mott for- 
mula (|6"6"|) and Wiedemann-Franz law (|68[) are not gen- 
erally applicable to systems where inelastic scattering of 
electrons plays an important role. Thus one cannot ex- 
pect to find the transport coefficients S, II, and K in 
long quantum wires by combining these relations with 
the expression ([60)) for the conductance. Below we find 
the transport coefficients of a long wire, whose length 
L ~ Z cq oc e^t 7 '. We will see that in such wires the 
thermopower and Peltier coefficient are no longer expo- 
nentially small, whereas the thermal conductance K is 
suppressed at L — ► oo. 

Unlike the thermopower and thermal conductance, the 
Peltier coefficient IT is defined in a system to which no 
temperature bias is applied, see Eq. (|63p . One can there- 
fore obtain II using the results of Sec. IIIII In particular, 
we saw that in a long wire the heat current is deter- 
mined by the parameter u of the distribution function, 
see Eq. (|48|). The value of u depends on the length of 
the wire via expressions ([57| and (|58|) . Combining these 
results we find the heat current in the form 



7T 2 T 



L 



6 fx l cq + L 



nv d . 



The ratio of jq and the electric current / 
the Peltier coefficient 



(69) 



envd gives 



6e fi l c 



(70) 



Similar to our main result (|60[) for the conductance, 
Eq. (I70|) is applicable at L l\. It shows how II grows 
from exponentially small values at L <C l C q to 7r 2 T 2 /6e/i 
at L —* oo. 



The thermopower and thermal conductance are defined 
in a system with a small temperature bias AT. To find 
these transport coefficients we revise our analysis of the 
conservation laws (l4"Tj) and (j¥5|) to add finite AT. The 
left-hand side of Eq. (|4Tj) represents the particle current 
j = I/e in the wire with electron distribution (I61| . As 
we saw above, the thermopower of such a wire (|65p is 
exponentially small, and thus the effect of temperature 
bias on the current j is negligible. We thus have j = 
2eV/h and recover Eq. igg) . 

The left-hand side of Eq. (j45]) is the heat current in 
a wire with electron distribution (|61[) . which does not 
vanish at AT ^ 0. It is determined by the thermal con- 
ductance (I67[) of a non-interacting wire. Thus instead of 
Eq. (g7|) we obtain 



27T 2 

3h 



■TAT 



JQ 



(71) 



The right-hand sides of equations and ([7T|) are 
not directly related to the voltage and temperature bias 
of the wire, but are determined by parameters A/x and 
u of the partially equilibrated distribution (|35[) . Indeed, 
the currents j and jq do not depend on position, and 
can be calculated for the internal region of the wire us- 
ing Eqs. (|3"6"|) and (|^5|) . The backscattering of the right- 
moving electrons predominantly happens inside the wire, 
at distances over l\ from the leads, where the partially 
equilibrated distribution ()35|1 is established. Thus we can 
express N R and Q R in terms of A/i using the results |55|) 
and (|5D|) of Sec. IIIII As a result, we obtain the following 
two linear equations upon the parameters Afj, and u, 



2A/i 

(1 + fx)— h nu = 



-n- 



2An 



TqUU 



2eV 
~' 

7T 2 T 

3/i /i 



-AT. 



(72a) 
(72b) 



The system of equations (|72p can be easily solved. Then, 
substituting the resulting Afi(V, AT) and u(V, AT) into 
Eqs. (|3"6"|) and ([^5]) . one finds the electric current I and 
heat current j'q. 

On the other hand, it is easier to find the thermopower 
S and thermal conductance K by noticing that their defi- 
nitions (|62"|) and ([64]) assume the condition of zero current 
/. Then, from Eq. (|36|) we immediately find = — nu 
and equations ([72|) reduce to 



— r\nu = 
(ri + r )nu — 



2eV 

7T 2 T 

3/i \i 



-AT. 



(73a) 
(73b) 



Excluding the unknown parameter u, we find the linear 
relation (f6"2"|) between V and AT with 



S 



7T 2 T L 



6e /i Z cq + L 



(74) 
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Predictably, the thermopower (|74|) and the Peltier coef- 
ficient ([70]) satisfy the Onsager relation II = ST. 

Furthermore, using Eq. (|73bj) we express the heat cur- 
rent PS]) in terms of AT as 



2n 2 T 



JQ 



ro 



-AT. 



3h n + r 

Thus the thermal conductance takes the form 
v 2ir 2 T l cq 



3h l cq + L 



(75) 



(76) 



At L <C l C q, Eq. (176")) recovers the result (f6"7| for nonin- 
teracting wires, but as the length of the wire grows, K is 
suppressed as 1/L. 

The fact that the thermal conductance K vanishes at 
L — > oo can be understood as follows. In an infinitely 
long wire the distribution function of electrons reaches 
the fully equilibrated form ([3|) controlled by three pa- 
rameters, T, n, and the drift velocity u<j. The thermal 
conductance is defined under the condition that the elec- 
tric current / = envd vanishes, see Eq. (|64|) . Thus va = 
and the distribution ^ takes the form of the standard 
Fermi-Dirac distribution. Due to its symmetry p — > — p, 
the heat current jq vanishes, regardless of the temper- 
ature bias AT applied to the wire. One therefore finds 
that in an infinitely long wire K = 0. 

The thermoelectric properties of a device are some- 
times summarized in the form of the dimensionless figure 
of merit defined as 



ZT ' = 



GS 2 T 



K 



(77) 



The figure of merit measures the efficiency of thermo- 
electric refrigerators. As ZT diverges, the device attains 
Carnot efficiency. For a material to be a good thermoelec- 
tric cooler, it must have a high value for ZT and a typi- 
cal figure of merit ZT ~ 3 would make solid-state home 
refrigerators economically competitive with compressor- 
based refrigerators. 41 However, in many materials the fig- 
ure of merit is limited by the Wiedemann-Franz law, and 
remains near 1. 

Substituting our results (|74|) and (|76|) for a long wire 
and neglecting a small correction to the quantized con- 
ductance, we find 



TT 2 (T^ 

" 12 \fij l eq (l e 



(78) 



In short wires the thermopower is small, resulting in a 
small ZT . On the other hand, the thermal conductance 
K is strongly suppressed in long wires, giving rise to in- 
finite figure of merit at L — > 00. 



V. DISCUSSION 

In this paper we studied the transport properties of a 
partially equilibrated quantum wire. In one-dimensional 
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FIG. 6: Correction to the conductance of a quantum wire as 
a function of its length. For wires longer than the scale h 
the conductance is given by Eq. (|60[1 . The behavior of the 
conductance in the regions L <C Iq and Iq <C L l\ is given 
by Eqs. (|6} and (|25[) . respectively. 



systems, equilibration of weakly interacting electrons is 
strongly suppressed at low temperatures, and the re- 
sulting equilibration length Z oq is exponentially large, 
Eq. (|59|) . Our main result is the expression (|60p for 
the conductance of a wire whose length L exceeds the 
length scale l\ given by Eq. (|32p . Because the scale l\ is 
only power-law large at low temperature, the expression 
([60)1 describes the full crossover behavior of conductance 
between the regimes of negligible and full equilibration. 
We have also been able to establish a connection between 
our result (|60p and the expression ^ for the correction 
to the conductance of a short wire obtained by Lunde 
et aZ3I. Similar to Eq. ©, our result (|6"U|) is exponen- 
tially suppressed at small L and grows linearly with L. 
However the prefactors are parametrically different. This 
mismatch is resolved by noticing that Eq. ([6]) is valid 
at L <C lo, where the length Iq defined by Eq. (|14p is 
short compared to l\. In the regime of intermediate wire 
lengths, l <C L <C h, the correction to the conductance 
(|25f scales with the length as L?l z . A summary of our re- 
sults for the conductance of a quantum wire as a function 
of its length is presented in Fig. 

In addition to conductance, we studied thermoelectric 
effects in partially equilibrated wires, limiting ourselves 
to the most interesting regime L ^> l±. The equilibration 
of the electron system has a dramatic effect on the ther- 
mopower and thermal conductance of the wire. As the 
length of the wire increases, the thermopower increases 
dramatically, from exponentially small values at L <C ^ C q 
to S ~ T/efi at L >• l eq , see Eq. ([73]). Conversely, the 
thermal conductance of the wire decreases due to the 
equilibration of the electron system in the wire from the 
Wiedemann-Franz value K — 2n 2 T/3h to zero, Eq. (|76|) . 



As a result, at L > ^ oq the quantum wire becomes a 
perfect thermoelectric refrigerator . 

In this paper we accounted for the effect of electron- 
electron interactions but neglected the electron-phonon 
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scattering, which may also affect the electron distribu- 
tion function. In GaAs quantum wires the phonons are 
three-dimensional and in equilibrium with the rest of the 
system. Scattering of electrons by phonons should there- 
fore have the effect of equilibrating them in the station- 
ary reference frame, thereby reducing the degree of equi- 
libration ?/. We leave the detailed study of the effect 
of phonon coupling to future work and limit our dis- 
cussion here to a few qualitative remarks. The reason 
the electron-phonon coupling is typically neglected com- 
pared to electron-electron interactions is that respective 
coupling constant is much smaller for the phonons. On 
the other hand, we do not expect the effect of phonons 
on the electron distribution function to be exponentially 
suppressed as e~^/ T . Thus we expect that coupling to 
the phonons to be negligible only at not too low temper- 
atures. 

Another effect neglected in this paper is the possible 
presence of slight long-range inhomogeneities in the wire, 
which would typically be caused by the presence of re- 
mote impurities in the GaAs heterostructure. The effect 
of such inhomogeneities on the conductance was studied 
earlier— 12£ under the assumption of full equilibration of 
the electron system. On the other hand, the inhomo- 
geneities themselves resist the equilibration process, and 
we expect an interesting interplay of these effects in the 
case of a partially equilibrated wire. We leave such a 
study to future work. 
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Then the homogeneous component Vq of the interaction 
potential takes the form 



(A2) 



For small wave vectors q <C 1/d, the Fourier trans- 
formed potential V q departs from the homogeneous com- 
ponent Vq by an amount 



V Q -V 



dx [cos(<7x) 
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2e 2 
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It then follows that the small-g behavior of the Fourier- 
transformed potential is given by 



V q = V 



1 - q z d 



2 MU\Q\d) 



\n(d/w) 



(A4) 



This expression contains an extra logarithmic-in-g factor 
compared to the expression introduced by Lunde et al.^- 



V„ = Vq 



(A5) 



However, as argued in the text, the typical scattering 
processes studied here only involve small momentum ex- 
changes, of the order hq ~ T/vf, so that the expression 
4j for V q reduces to the one of Eq. (|A5|) with 
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The model introduced in Eq. (|A1[) . therefore merely 
amounts to an extra logarithmic temperature dependence 
in the length l eee , Eq. (TTTj) . 



APPENDIX B: FOKKER-PLANCK EQUATION 
FOR 3-PARTICLE COLLISIONS 



APPENDIX A: SCREENED COULOMB 
INTERACTION 

Let us consider the Coulomb interaction between elec- 
trons screened by a nearby gate, which we model by a 
conducting plane at a distance d from the wire. In this 
case, the electron-electron interaction takes the form 



In this appendix we discuss the Fokker-Planck approx- 
imation and calculate the coefficients for the interaction 
potential used by Lunde et alM- as well as the screened 
and unscreened Coulomb interaction. 



1. Fokker-Planck approximation 



V(x) 



1 



1 



M ^/x 2 + (2d) 2 



(Al) 



The diverging short-range behavior of this potential 
needs to be regularized in order to evaluate the small- 
momentum Fourier components V q . To this end, we in- 
troduce the small width w of the quantum wire, w <C d. 



We start out from the collision integral © for the 
three-particle scattering process. We discussed in sec- 
tion |TT] that the only contributions relevant to transport 
result from collisions involving two pairs of incoming and 
outgoing states in the vicinity of the right and left Fermi- 
points, and one pair of incoming and outgoing states at 
the bottom of the band, see Fig. [2Jd. Let pi and p[ be 
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the momenta near the bottom of the band, p 2 and p' 2 
the ones near the left Fermi point, while p 3 and p' 3 are 
taken near the right Fermi point. Unprimed momenta 
correspond to incoming states whereas primed ones are 
associated with outgoing states. We introduce the hole 
distribution g Pi = 1 — f Pi and the collision integral of 
holes I p ,x[g] = —Ip,x[f], which using Eq. (J5J) , can be re- 
cast as 



W[<?] = E [W(pi,p'i)9p{ - W(p' 1 ,p 1 )g Pl ] , (Bl) 



where 



W(pi,pi) = 48 J2 Ww.vvygvgvfxhh, (B2) 

P2,P3 
P2,P3 

Wfa,pi) = 48 W i3A ,v2>v92toh>fvfo. (B3) 



P2,P3 



W(pi,Pi) is the rate for a transition in which a hole scat- 
ters from some state p\ into px, while W(pi,pi) denotes 
the corresponding transition rate for the inverse process. 
Here and in what follows, all momentum summations are 
restricted to the ranges discussed above. This restric- 
tion results in a combinatorial factor of 12 in Eqs. (jB2|) 
and (|B3p . The remaining factor of 4 originates from the 
spin summations as we anticipated that the main con- 
tribution to the 3-particle scattering rate u>i23 ; i'2'3' of 
Eq. ijHJl takes the form 5 (Tl(7l ,S (T2rT2 ,8 (73(Ts ,Wx23;X'2'3' , with 
a spin-independent Wx23;V2'3' ■ This simplification is only 
valid in the limit of small momentum exchanges, and 
can be performed here since for the Coulomb interaction 
Vo ^> Vk F ■ Since px and p\ lie near the bottom of the 
band, the distribution functions g Pl and g p t are expo- 
nentially small, and so is the collision integral of holes 
(|Blj) . It is therefore unnecessary to account for addi- 
tional exponentially small contributions in the scattering 
rates W(px,p[) and W(p' 1 ,px), so that one can safely re- 
place /1 ~ 1 and f v ~ 1 in Eqs. (jB2j> and (|B3]) . 

The Fokker-Planck approximation exploits the fact 
that collisions typically induce small momentum changes 
of order 0{T/vp). For the following, it is convenient 
to introduce the momentum exchanges qi = p[ — p\ . 
With this notation, W(p[,px) describes the transition 
rate for the process in which a hole scatters with mo- 
mentum transfer qi, from the initial state jpi, and can 
thus be rewritten as W(p[,pi) = W qi (pi). Following the 
same prescription, the transition rate for the inverse pro- 
cess becomes W(pi,p' 1 ) = W- qi (pi + qi). Performing a 
small-momentum expansion, one has 



It 



Oi + gi)fl P 'i = w -n (pi)9pi + qid P1 (W- qi (pi)g pi ) 

+ f d pi (W-n (px)9 Pl ) + O {qld 3 pi (Wg)) , (B4) 



where d Pi — d/dpi. Introducing further 

91 91 
9l 91 

the collision integral of holes takes the simplified form 

hxM = ~ d Pl (A(pi)g(pi)) + \d 2 pi (B(px)g(px)) ■ 

(B7) 

We next turn to the explicit derivation of the functions 
A(p) and B(p) in the case of three-particle collisions. 



2. Relation between A(p) and B(p) 

The scattering rate Wi23;i'2'3' contains both the energy 
and momentum conservation and can be rewritten as 

Wi23;i'2'3' =8(ei> -ei + e 2 > - e 2 + e 3 > - £3) 

x ^91+92+93,0 w(qi,q 2 ,q 3 ), (B8) 

where we introduced the momentum transfers qi = p[ — 
Pi. The function w that remains after writing the con- 
servation laws explicitly, should depend on all pi and py . 
However, for the momentum configuration under consid- 
eration, p\ lies near the bottom of the band, while p 2 
and P3 lie near the left and right Fermi points, all within 
a small range set by temperature. We thus argue that, 
up to small corrections in T/ fx, one can replace pi ~ 0, 
p 2 — ~pf and P3 ~ pf in the expression for w, which 
then becomes a function of q\ , qi and q 3 . 

Using the approximated forms £2' — £2 — — vpq 2 and 
ey — £3 ~ vpq3, the conservation laws allow us to express 
q2 and q 3 in terms of px and p[ as 



92 = 
93 



Pi - Pi ty ~ ei 
2 2v F ' 

Pi - Pi ei> - ei 



2v F ' 



(B9) 
(BIO) 



where one readily sees that (72 — 93 — — 9i/2, up to small 
contributions of order px/pp <C 1. 

Substituting the expression (|B8[) for the scattering rate 
into Eq. (|B3|1 , and using the energy and momentum con- 
servation laws to simplify two of the momentum summa- 
tions, one has 

W(j)'i,px) = 48^w(<?i, q 2 ,q 3 ) ^ ^2^3/^2+92/^3+93: 

F P2,P3 

(Bll) 



where we focused on a section of the wire, of length AL <C 
l eq . Here q 2 and q 3 are functions of p\ and p[, as given 
by Eqs. JH§} and (|BT0| . 
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The remaining momentum summations can be per- 
formed explicitly upon linearizing the dispersion near the 
Fermi level 

X- , AL q 2 e^ 

9p2 Jp2+q-2 — u 1 „/)«,.«, ' \ olz ') 



P2 

$P3./p3 + 93 



hi- el 3vF ^ 

AL g 3 

hi- e^Fiz 



(B13) 



so that the transition rate, Eq. IjBllj) becomes 



(B14) 



V F V h J sinh 2 f§v^HL\ 



where we replaced q 2 and q 3 following Eqs. (|B9j) and 
(IB10[) . and introduced w qi = w(qx,—qi/2,—qi/2). The 
leading contribution to W(p'i,pi) is an even function of 
qi which leads to a vanishing A(pi) once substituted into 
Eq. (|B5|) . For that reason, we expanded the expression 
for the transition rate up to linear order in the small 
parameter (ey — £\)/T ~ ViIvf "C 1. 

The functions A(p\) and B(p\) are then readily ob- 
tained from Eqs. (IB5|) and (|B6|1 by substituting the ex- 
pression (|B14[) above for W(p'i,px) yielding 

A(pi) = 1^2<H \W{pi +qi,Pi) - W{pi - q%,px)] 



Pi 



3 / AL\ 3 qf w qi 



^7 I ~/T 
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and 



B (Pi) = \^2 < ll [W(pi + qi,Pi) +W(px - qi, Pl ) 



E 



3 / AL\ 3 qf w q 



V F \ h ) sinh 2 (^2L 



(B16) 



where we discarded contributions of order (j>i/pf) 2 and 
higher. It follows from these two expressions that for a 
momentum p\ deep in the band, \p\ \ <C pf, the function 
B(pi) can be approximated by a constant B, while A(p\) 
satisfies 



A(pi 



Pi 

2mT~ 



B. 



(B17) 



In order to derive an explicit form of the Fokker-Planck 
equation, it is therefore sufficient to calculate the con- 
stant B. 

Let us now briefly comment on the validity of the 
Fokker-Planck approximation. The first two terms ne- 
glected within the Fokker-Planck approximation would 
contribute to the collision integral as dp 1 (C(pi)g Pl ) 
and d^{D(px)g pl ), where C(pi) = J2 Ql Qi w qi(Pi) and 



D{pi) — J2 qi qf W qi (pi). Going through the same deriva- 
tion as the one outlined above, and keeping in mind 
that every new power of q results in a factor of T/vf, 

one can convince oneself that C(p\) oc f^rj MjP\) and 

D{pi) oc (^~\ B. It results that the contribution to 

the collision integral from the terms in C(p\) and -D(pi) 
are smaller than the ones from A(p\) and B(pi) by a 
factor T/fi <C 1. This readily generalizes to higher or- 
der derivatives <9™ , thus validating the expansion of the 
collision integral used here. 



3. Evaluation of B 

We now derive the expression for the constant B us- 
ing the specific form of the electron-electron interaction 
potential of Eq. (|A5[) . This expression of the potential is 
largest for small wave vectors q allowing to discard the 
exchange terms |Vjqi^fc F | <§C |Vj g |<gfc F | in the scattering 
rate, which is thus dominated by the direct term. Fol- 
lowing Ref . [27| , the reduced scattering rate w q takes the 
form 
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V q /2 (V q -V q /2) 

2AL 2 n 



(B18) 



Expanding for small values of q, this can be further 
rewritten as 



1 



_ 9tt2 

' Jq ~ 16/i (fc F AL) 4 ~ 



(B19) 



where Vq is the zero-momentum Fourier component of 
the potential. Substituting this expression back into 
Eq. (|B16[) , and performing the integral over qi , we finally 
find for the constant B 
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kpPFfJ; (B20) 



from which we can extract the length scales l\, l eee and 
l using Eqs.®, and (TTJ} respectively 



, fV k 



40 



qo 



i 15/2 



r 1 



In' 



Vokp\ 4 f kp\ A /T x 7 



qo 



(i 



5o 



k F , (B21) 
(B22) 
(B23) 



The expression for B is model-specific, and the result of 
Eq. (|B20p was obtained for the potential V q of Eq. (|A5[) . 
leading to the same expression for l eee as the one of Lundc 
et al.r^ Eq. (jlip . In the case of the screened Coulomb 
potential discussed in Appendix [X] and for temperatures 
T <C hvp/d, one obtains similar results upon redefining 
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go according to Eq. (|A6|) . On the other hand, for tem- 
peratures T ^> hvp/d, the effect of the screening gate 
can be neglected, and the electron-electron interaction is 
then well described by an unscreened Coulomb potential, 
of the form V q = — ln(l/\q\w) at small wave vectors 
\q\ <C This in turn leads to the following value of B 

(B24) 

and the corresponding expressions for the length scales 



, 4^Fln 2 2 fe 2 k F \ i /T 



e /J, 



7/2 



In 2 f ^ ) /,-,. 

1 w 



r 1 



e^kp \ ' IT 
e n 

e 2 k F 
e \x 



In 2 



hvp 

Tw 



kF, 



/,T^ (-^) (-) ln 2 (^)^- 



(B25) 
(B26) 

(B27) 



Substituting the expression (|A2[) for Vo, and (|A6[) for 
qo into Eq. (|B20|1 . one readily recovers that Eqs. (|B20[) - 
(|B23|) match with Eqs. (|B24|) - (|B27| at the crossover tem- 
perature T ~ hvp/d. 



APPENDIX C: POSITION DEPENDENCE OF 
THE DISTRIBUTION FUNCTION 

In this appendix we determine the profile of the 
position-dependent parameters /i R ^ L (x) and T{x) enter- 
ing the distribution function (|35|) . 

As discussed in the text, steady state parameters u 
and A/t are constant along the wire. Therefore, the spa- 
tial profile of distribution (|35[) is determined by space- 
dependencies of the average chemical potential and tem- 
perature. It is convenient to measure deviations from the 
lead values, 



T(x) = T + 5t{x), 



(CI) 



and /i L (x) = /i R (x) — A/i. Let us then consider a wire of 
length L>!i, and focus on a small segment between the 
positions x and x + Ax, where < x < L. We observe 
that conservation of momentum insures homogeneity of 
the momentum current 



Jp 



3 P + ^(eV + 6n R + 5 f i L ) + ~p F ST, (C2) 



A 4 



where j P is the momentum current in absence of external 
potential bias. From Eq. (|C2[) one readily observes that 
for jp to remain constant, a drop in chemical potentials 
must be compensated for by an increase in temperature, 



5t(x + Ax) — St(x) 
5[i R {x + Ax) - 8ii R {x) 



6 pL 



(C3) 



valid up to small corrections in (T//i) 2 . To calculate the 
spatial profile of /j, r and T we need to find the slope of 
either one and their boundary values near the leads. 

The slope of (jl r is readily found from calculating the 
difference in right-moving particle currents within the 
segment of length Ax, 

j R (x + Ax) - j R (x) = | {5n R {x + Ax) - 6v R (x)) . 

(C4) 

Since this difference equals the rate of change of the 
number of right-moving electrons within the segment 
(x,x + Ax), we can insert Eq. (|55|) into the left hand 
side, to find 



d/i R (x + Ax)-Sfi R (x) 



hh R 



A.i 



where n R = N R /L. Then, from Eq. ([55 



dfj, R {x) 
dx 



eV 



12 \n J L + l. 



(C5) 



(C6) 



eq 



i.e., the chemical potentials linearly decrease, while the 
temperature linearly increases along the wire, 



» L (x)= f i L (L) + 



eV x 
12 \/i J L + l cq ' 

eV(L - x) 



^2 /rp\ 2 



T{x) = T(0) + 



12 \n 
T eVx 
2/7 L + l 



L + l 



(C7) 

(C8) 
(C9) 



eq 



In the linear response regime boundary values n R (0) 
and H L {L) deviate from the chemical potentials in the 
leads by an amount proportional to the applied voltage. 
From inversion symmetry it further follows that these 
deviations are opposite in sign, i.e. 

M R (0) =fH - XeV, (CIO) 
H L {L) =fi r + XeV. (Cll) 

The parameter A may be inferred from the equation 

fi R (Q) - ji L {L) - fi R (0) - fi R (L) - A/x, (C12) 



by inserting Eqs. (|C10[) and (jClip into the left hand side 
and rewriting the right hand side using Eqs. (|C6[) . (|5T 
and (jnnj) ■ This results in 



A = 



hnu 



(C13) 



where n is the electron density. 

The boundary values T(0) and T(L), are found in a 
similar way by combining Eq. (|C9|) with the observation 
that (5r(0) = — 8t{L), which is, again, a consequence 
of inversion symmetry. We summarize the values for the 
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position-dependent parameters close to the leads in terms 
of eV and u 



H R {0) = fi + eV - 



hnu 



V L {L)=n + 



hnu 



T(0)=T(1--), 



T(L) = T 1 + — , (C14) 



where vp — -\/2/x/m. 

Finally, restricting ourselves to linear terms in V and 
finite-temperature corrections to leading order in (T//1) 2 , 
we observe that the values given in (|CJ14|) guarantee that 
all moments (v p p s ), with s € N, are continuous at the 
boundary between the wire and the leads, i.e. 



dp 



v p p 



dp 



/ P (o)-4 0) 



f P (L)-6 0) 



= o, 

= 0, 



(C15) 
(C16) 



where fp ^ is the lead distribution function |T]). For s = 0, 
1 and 2, relations (|C15|) and (|C16|) imply continuity of 
particle, momentum and energy currents at the boundary 
between wire and leads. 

To show the validity of Eqs. (jCl~5)) and (|UT6|) we ex- 
press the distribution function (|3"5)) in terms of 8[i RL and 
8t, and expand the difference of distributions entering 
Eqs. (|C15p and (|C16p to linear order in these parame- 
ters. For right-moving electrons close to the left lead, 
one has 



/ P (o)-4 0) = - 



5/j, r (0) + up+ (e p - fi) 



St(0) 



T 



(0) 



de 1 
(C17) 



and similarly for left-movers close to the right lead. Upon 
introducing the new variables £ = p 2 /2m — /i and z — 
£/T, and neglecting exponentially small contributions oc 
g-W T thjg results in 



dp 

T v pP 



ip(0)-4 o) 



(2m/x.): 



zT 

dz [ H 



(5^(0) + u^2m f Jl + ^ + St (0)z 



M, (C18) 



where fa(z) = (1 + e z ) . Keeping now only terms up to quadratic order in (T/fi) one then finds 



dp 



v p p 



/ P (o)-/i 0) 



(2m,fj,)'- 



Sf, R (0) 



2rp2l 



l + s(s-2) 



n 2 T 
24 /i 2 



+ uy2m/i 



2 ^21 



l + (.s 2 -l) 



St(0)s- 



n 2 T 
6 /i 
(C19) 



Substituting values for (5/i fl (0) and (5r(0) from (|C14[) . one 
can readily check that Eq. (|C15|) is satisfied for all values 
of s E N. Proceeding in an analogous way at the right 
end of the wire confirms Eq. (|C16|) . 



APPENDIX D: ENERGY TRANSFERRED IN A 
BACKSCATTERING PROCESS 



In this appendix we calculate the change in right- 
movers energy associated with the backscattering of a 
right-moving electron. 

Let us focus on a small segment of wire in between 
positions x and x + Ax. Following Eq. (f39|) . we use the 
conservation of the number of particles to express the rate 
of change in the number of right-movers N R in terms of 
particle currents. Proceeding similarly with the energy 



currents, one can express the ratio E R / N R as 

E*_ = j R (x + Ax) ~ j R (x) 
NR j R {x + Ax) ~ j R (x) 

ir 2 T 5t(x + Ax) - St(x) 
3 &ii R {x + Ax)~5iJL R {xy ( } 

where we used the distribution function (f35|) to calcu- 
late the current differences in terms of St and S[i R . The 
first contribution fj, is the energy carried by the elec- 
tron making its transition from the subsystem of right- 
to that of left-movers. The second contribution repre- 
sents the energy of excitations created at the right Fermi 
point during the sequence of three-particle scattering pro- 
cesses that ultimately results in the backscattering of a 
right-mover. This contribution can also be viewed as 
the heat transferred from the right- moving subsystem for 
each backscattering process, in which case, the prefactor 
ir 2 T/3 is readily obtained from the thermal conductance 
([67)) and the relation (|C4[) between right-movers current 
and chemical potential. 

Substituting the ratio of changes in temperature and 
chemical potential as given by Eq. (|C3[) into Eq. (|D1[) . 
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one has 



E R 
~NR 



1 + 



(D2) 



where the term in O (T/fi) is discarded in the text, as it 
only leads to subleading corrections to the conductance. 
We also briefly mention that Eq. (|D2[) was derived for a 



quadratic dispersion, but can be generalized to the case 
e P °c H s , yielding 



(1 - s)/x, 



(D3) 



again, up to subleading corrections in T / fi. 
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